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We consider the eigenvalue problem for one- dimensional linear Schrodinger lattices (tight- 
binding) with an embedded few-sites linear or nonlinear, Hamiltonian or non- conservative defect 
(an oligomer). Such a problem arises when considering scattering states in the presence of (gen- 
erally complex) impurities as well as in the stability analysis of nonlinear waves. We describe a 
general approach based on a matching of solutions of the linear portions of the lattice at the loca- 
tion of the oligomer defect. As specific examples we discuss both linear and nonlinear, Hamiltonian 
and VT- symmetric dimers and trimers. In the linear case, this approach provides us a handle for 
semi-analytically computing the spectrum [this amounts to the solution of a polynomial equation]. 
In the nonlinear case, it enables the computation of the linearization spectrum around the station- 
ary solutions. The calculations showcase the oscillatory instabilities that strongly nonlinear states 
typically manifest. 

We consider the time evolution of a quantum mechanical wave function as governed by 
the Schrodinger equation. The wave function is distributed spatially on a discrete one- 
dimensional lattice, i.e. a chain of nodes indexed by integers, so that the spatial derivatives 
are replaced by differences. The potential function is nonzero only at a few center sites on 
the lattice, representing either physical impurities or other obstacles such as an external 
field or a nonlinear material. Since the general solution of the zero potential problem is 
well-known, we begin by constructing it on the outer (left and right) portions of the lattice. 
Working our way toward the impurity sites using the restraints of the discrete Schrodinger 
equation, we find that appropriately defined portions of the outer solution must satisfy a 
polynomial equation. Our method is also applied to the (again) discrete but (now) nonlinear 
Schrodinger equation. Here known stationary solutions are acted upon by a time-dependent 
perturbation, and we find that appropriately defined portions of the perturbation must 
satisfy polynomial equations. The main point is to show that the polynomial conditions 
we derive accurately determine the dynamical stability of the solutions. The success of our 
method in tracking the associated linear and nonlinear spectra is presented throughout the 
linear and nonlinear cases. In order to demonstrate the generality of our approach we show 
examples using both real valued Hamiltonians and purely complex parity-time symmetric 
potentials. 

I. INTRODUCTION 



Linear lattices with embedded impurities have been considered to model many different physical systems. A 
familiar example is the case of mass defects in an otherwise pure harmonic crystal. In the context of the tight- 
binding description of electron transport, the nonlinear terms describe the strong interaction with local vibrations at 
the impurity site Models of the same type have been used to describe tunneling through a magnetic impurity 
connected to two perfect leads in the presence of a magnetic field [2]. Another wide domain of application is the 
one of nonlinear optics. The type of system thereby considered is an array of coupled waveguides consisting of (at 
least) two types of materials one linear and one of Kerr type, see e.g. [3|. Several problems have been addressed, 
ranging from the existence of localized solutions residing on the nonlinear portion (nonlinear impurity modes [4|) to 
the scattering and transmission of plane waves through it |5j. 

An important question concerns the spectral characteristics of the problem in the linear regime and the dynamical 
stability of the solutions in the nonlinear regime (i.e., when the embedded impurities are genuinely nonlinear). To 
the best of our knowledge no systematic study of such an important issue been presented in the literature. In Ref. [(| 
the continuum case of Nonlinear Schrodinger equation with a localized (^-function) nonlinearity was considered. 
For discrete lattices, an analysis of a related problem, the nonlinear Fano effect, has been reported in [7] (see also 
references therein). Another related reference is [8] where a bifurcation analysis of an open chain with nonlinearity 
and disorder is performed. 

The aim of the present paper is to formulate (in general) and to solve (in some specific cases) the spectral problem 



n 2 

for a discrete Schrodinger equation with TV embedded defects (an TV-site oligomer [9[). Equations of this type emerge 
in, at least, two interesting cases. The first is the one of linear chains with localized complex on-site potentials. The 
second one arises when considering the stability problem of complex (i.e. current carrying) scattering solutions of the 
discrete nonlinear Schrodinger equation [E| . Such solutions can be computed exactly by the transfer matrix method 
||J fTol - H^I and linearization around them yields linear equations of the same type. 

Our motivation to approach this issue is twofold. The first is to study the scattering problem by open systems 
possessing a T^T-symmetry [l3l-[l5j. This theme is rapidly evolving into a major area of research partly due to its 
providing an intriguing alternative set of non-Hermitian Hamiltonians with potentially rea l eigenv alues and partly 
due to its experimental realizations in the field of nonlinear optics, both in the case of dimers [16|,[l7j and even in that 
of whole lattices [18| . The second motivation stems from a recent study of non-reciprocal wave propagation through 
non-mirror symmetric nonlinear lattices [19|. These are a minimal model for a "wave diode", namely for devices 
capable to selectively rectify wave energy through a nonlinear medium of propagation. Relevant applications arise 
for instance for phonon scattering at a nonlinear interface layer between two very dissimilar crystals [2pj |. acoustic 
waves in sonic materials [21] or in the so-called all-optical diode for photonic applications (22j . 

The solution of infinite-dimensional problems like the above is technically more involved than a straightforward 
matrix diagonalization. Indeed, to solve the problem exactly one should impose the solution a definite plane wave 
form (with generally complex wave numbers) in the two semi-infinite linear parts of the chain. The matching of such 
waves through the oligomer portion reduces the infinite-dimensional problem to an homogeneous linear system of 
2TV equations, whose solvability condition, along with the dispersion relations, yields a set of nonlinear equations for 
the unknowns. 

The details of this method will be presented for a linear PT-symmetric case [l3l-[l5|. In the latter, the imposition 
of the parity (V, associated with spatial reflection) and time-reversal (T, associated with temporal reflection and also 
i —> —i) symmetries leads to an imaginary (growth/decay) part of the Schrodinger potential which is spatially odd. 
The relevant more general constraint, in fact, reads V(x) = V*(— x) (where ★ stands for complex conjugation); i.e., 
the real part of the relevant potential needs to be even. This setting will be compared/contrasted to the Hamiltonian 
case, where the potential is real. 

The paper is organized as follows. In section [III we will discuss the simplest case namely that of linear oligomers. 
Although simple, this case is useful to set the stage of the general methodology and to illustrate some basic features 
that will be useful throughout the rest of the paper. In section ITTTl we turn to the nonlinear case. In each of these, 
we first present the relevant analysis and subsequently we corroborate it on the basis of numerical computations. 
Finally, in section [TV] we present a summary of our findings and some potential directions for future study. 



II. LINEAR CASE 



In this section we present the strategy to compute the spectrum and eigenstates of a general linear problem where 
TV-site complex defects are embedded in an otherwise homogeneous lattice. More precisely, we consider the dynamic 
evolution of a one-dimensional chain governed by the linear discrete Schrodinger equation 

l(j) n - V n (j> n + n +l + (j> n -l = (1) 

where 4> n (t) G C. The complex potential V n G C is zero everywhere except for 1 < n < TV for some integer TV > 2. 
We look for solutions of the form <j) n (t) = a n e %UJt so that (pQ) becomes the time-independent condition 

Gn+l + - (V n + u)a n 0. (2) 

Note that a n is complex as well. 



A. Theoretical Analysis 

Our goal is to present a strategy for the analytic computation of the eigenstates a n and eigenvalues lj, assuming 
that V n has been prescribed. First we begin with an observation which applies to the portions of the lattice on the 
left and right of the TV central sites. The zero potential (V n = 0) version of equation ([2]) has a general solution of the 
form a n = Ae lKn + Be~ %Kn for A, B G C and k G C satisfying the dispersion relation 2cos(ft) = uo. Thus we begin 
with the ansatz 

an ~ \ Ce iKn + De~ iKn n> TV W 

for A, C, 5, D G C and k satisfying 



1 

uj = z H — 

z 



(4) 



where we use the shorthand notation z = e %K ^ 0. Imposing (j3]) and (j4]) is enough to satisfy (j2j) for all n except 
n = 0, 1, . . . , A", N + 1. Applying ([3]) to ([2]) with n = 0, A 7 " + 1 shows that the formula ([3]) is applicable at n = 1, A\ 
Since in practice we deal with a finite length lattice we introduce indexing no < n < mo for large magnitude integers 
no << 1 and tuq » N. Imposing [27] homogeneous Dirichlet boundary conditions a no _i = a mo+ i = shows that 
v = a no+ i/a no = a mo _i/a mo . Combining this with © gives B = -Ae 2i < nQ ~^ and D = -Ce 2i ^ mQ+1 \ Thus we 
have 

_ f Ac n n < 1 ( v 

for C n = Z n - z 2 ™o-2-n and ^ = _ ^mo+2-n_ 

At this point the imposition of (|5]) and ([!]) is enough to satisfy all but the N equations in ([2]) with n = 1, . . . , N. 
In the following sections we outline strategies for completing the computation of a n and uo in the cases of embedded 
few-site defects and apply the method to our P7~-symmetric or Hamiltonian oligomer. 



1. Oligomer of length two 



For N = 2 we examine the two equations ([2]) for n = 1 , 2 which by ([5]) can be written as 
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ai 
a 2 
a 3 
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C 



for P 



1 -(Vi+a;) 1 
1 -(F 2 +cj) 1 
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ci 

d 2 

^ 3 



(6) 



In other words, we now have a 2 x 2 system with unknowns A, C and coefficient matrix PQ. Imposing zero determinant 
det(PQ) = gives a solvability condition, which by ([!]) amounts to requiring that allowable z are roots of a polynomial 
with coefficients in terms of V n . For each root z one computes the corresponding A,C by the equation in (j6]) so 
that finally an eigenvalue uo and eigenvector a n are known by (|2|) and ([5]). Note that roots z = ±1 correspond to 
= a n = and should be disregarded so as to obtain non-trivial eigenvectors. 
As an example, consider the special case V\ + V2 = 0. One can compute by hand that the determinant condition, 
when expressed in terms of a positive power polynomial, is 

= ^2(mo-n )+6 + _ 1 ^2(m -n )+4 _ ^^2m +l _ y^y^mo 

_y 2Z 2m -l _ y 2Z -2n +7 _ ^^^-2^+6 _ y^-2n +5 + _ + j ( ? ) 

for which z = ±1 are each double roots. This reduces the number of relevant roots of Q to 2 (mo — no + 1) = 2L, 
twice the length of the lattice. One can easily show, using the symmetry of the coefficients in ([7j), that solutions of 
occur in reciprocal pairs z,K Combining this with (J2J) shows that our analytic computation indeed yields exactly 
L eigenvalues cj, counted with multiplicity. 



2. Oligomer of length three or greater 



For N > 2 the strategy is similar to that of the last section but involves an additional step. First focus on the 
N — 2 equations (|2j) for n = 2, . . . , N — 1 which we write as 



1 — (V2 + oj) 1 



1 







ai 


1 




a N 



0. 



-(Vjv-i+w) 

By moving ai in the first equation and a at in the last equation to the right-hand-side, J5J) can be written as a 
(A" — 2) x (A" — 2) system with unknowns . . . , a/v-i that can be computed linearly in terms of ai,a/v 
appropriate conditions of invertibility of the resulting coefficient matrix). For the case of N = 3, Eq. (j8j) 
equation which we solve for a 2 in terms of ai, (23 to obtain 

|5j c\A + d 3 C 



a 2 



for 7^ — V2. 

Continuing with the A" 



F 2 + V2 

3 case, we use (|5]) and (|9]) to rewrite the two remaining equations in ([2]) with n = 
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FIG. 1: The left set of panels corresponds to the linear VT- symmetric dimer with N = 2, ai,2 = and Vi = 27 = — V2 while 
the right set of panels to the case of the trimer with N = 3, 0:1,2,3 = and V± = i*y = — V3, V2 = 0. Each set contains a (top) 
plot of linear stability eigenfrequencies as a function of increasing 7 as computed an a lattice of length 200. The eigenvalue 
agreement between numeric (circles) and exact (x's) computations is shown for 7 = .5 (bottom left) and 7 = 1.5 (bottom 
right). 
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FIG. 2: Profiles of modulus of unstable eigenvectors \a n \ are shown for the linear dimer (left) and the linear trimer (right). 
These eigenvectors correspond to the eigenvalues with negative imaginary part for 7 = 1.5 which are seen in the bottom right 
panels of Fig. [1] Again, numerical computations (circles) are shown to agree with the exact computation (x's) and the lattice 
length is 200. 



Again the solvability condition is det(PQ) = so that z is a root of a polynomial with coefficients in terms of V n . 
For each z, one computes A, C in the nullspace of PQ from (fTUj) so that by (|3J), (|SJ) an d ©, an eigenvalue uj and the 
corresponding eigenvector a n are obtained. The N > 3 case is similar in that P and Q are determined using (J5j) and 
the above described process of using J5J) to obtain expressions for g^, ... , a/v-i m terms of A, C. 

Due to the condition obtained in (j9j), if there exists an eigenvalue uj = — V2 then it will not be found by the above 
process. In the case that such an eigenvalue exists the corresponding eigenvector is found by using (|5]) to rewrite the 
three n = 1,2,3 equations in (|2]) as the 2x2 system 

' (V 1 -V 2 )c 1 -c (V 2 -V 3 )d 3 + d 4 
ci d 3 

for z such that z + ^ = — V2. Then a n is computed by ([5]) and a 2 = (V\ — V 2 )ai — ao = (V3 — V 2 )a 3 — 0,4. 
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FIG. 3: The left set of plots corresponds to the linear Hamiltonian dimer with N = 2, ai,2 = and Vi = 1 = V2 while the 
right set of plots to the case of the linear trimer with N = 3, 0:1,2,3 = and V\ = 1 = V3, Vi = 0. Each set contains a (top) 
plot of linear stability eigenvalues and a (bottom) plot of the eigenvector for the defect (point spectrum) mode. All plots show 
agreement between the numerical (circles) and semi-analytical (x's) results. Here the lattice length is 40. 



B. Numerical results 



First, we focus on the PT-symmetric case where the dimer has linear potential V\ = — V2 = ij and the trimer 
potential is V\ = —V3 = 27, while V2 = 0. Fig. [1] shows agreement between the eigenvalues uo computed numerically 
directly from (|2]), through an eigenvalue solver, as compared to the semi-analytic calculation involving the identifi- 
cation of the roots of an equation of the form of Eq. ([7]), as outlined in the above theoretical analysis section. Fig. 
[2] shows the accuracy of the semi- analytic calculation in predicting the modulus profile of the eigenvectors \a n \. 

The spectra in Fig. [I] have both, as expected, a continuum component filling densely the interval [—2^2] on the 
real axis corresponding to the propagation of linear waves. For 7 > j c the solutions become unstable [28|; this is 
the so-called PT-phase transition of [l3l-[l5j. In the case of the dimer, the principal unstable eigenfrequency pair 
has zero real part and is purely imaginary. In the case of the trimer, on the other hand, there is some oscillatory 
behaviour superimposed to the exponential growth of perturbations. The exponential localization of the associated 
eigenvectors is associated with the localized nature of the embedded defect structure and is showcased in Fig. [2j 

Finally, as regards the linear eigenvalue problem, we consider the Hamiltonian case where the dimer has a linear 
potential Vi = V (l + 6), V 2 = V (l - S) and the trimer potential is V x = V (l + 6), V 2 = 0, V 3 = V (1-S). Fig. [3] 
shows agreement between the numerically computed and semi-analytically calculated eigenvalues and the (modulus 
of the localized) eigenvectors. A fundamental difference here concerns the Hermitian nature of the relevant (matrix) 
operator which excludes the possibility of imaginary eigenfrequencies. Nevertheless, there exists within the spectrum 
a real defect eigenvalue. For S > 0, this frequency decreases as S increases and the corresponding eigenvector profiles 
for 5 > are typically similar to what is seen in Fig. [3j 

III. NONLINEAR CASE 

Let us now consider the nonlinear Schrodinger equation [23| 

i$ n - V n $ n + $ n+ i + $ n -i = Oi n \^ n \ 2 ^ n (12) 

for a n G R zero everywhere except 1 < n < N. Seeking, as is customary, stationary solutions of the form <E> n (t) = 
^ n e~ %UJt : , we find that they should satisfy 

~ V n ^ n + + = OL n \lp n \ 2 1p n . (13) 

We consider the dynamics of small perturbations defined by 

*n(t) = OAn + 60 n )e-^; <t> n = (a n j vt + b n e~ %v * £ ) (14) 

for uo G IR and a n ,b ni v G C and with ifj n G C. In order to investigate the stability, we examine the resulting order-e 
equations, amounting to the spectral or linear stability analysis equations 

i<j>n = (V n - Uj)(j) n - n+ i - n _i + a n (2\^ n \ 2 (j) n + IpnK) ■ 



Note that (j) n is complex. It is thus recognized that problem ([T5]) is similar to (pQ), the main difference being that now 
<j) n is coupled to 0*. Hence, given a zeroth-order solution ip n , we can apply a similar approach as the one developed 
in section [III in order to determine its linearization spectrum. 

For later reference, it is convenient to reformulate the problem in matrix form for a n , b n which are obtained from 



12]) by equating coefficients of e 



i(y—u))t p — i(v* -\-uo)t 



yielding 
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F3 F4 



and 



Fl 
^3 



diag(w -V n - 2a n \ip n \ 2 ) 
diag(a„(C) 2 ) 



a n 
b* 
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F 2 
F 4 



a n 
b* 



-diag(a„^£) 
-diag(w - V* 



(16) 



2aM 2 )-G 



(17) 



where G is a sparse matrix with ones on the first super- and sub-diagonal. 



A. Theoretical Analysis 

Similar to the linear case we begin with the linearization problem ansatz 

r Ae %Kn + Be -*«n n<1 C A'e iK ' n + B'e- iK ' n n< 1 ( . 

an ~ \ Ce iKn + De~ iKn n> N ' n ~ \ C"e™' n + D'e^' 71 n> N ^ ' 

for A, C, 5, D, A', C", P', D' G C and k, G C satisfying dispersion relations 

z^ — cj = 2cos(^) — [y + cj) = 2 cos(V). (19) 

Imposing ([18]) and ([T9]) is enough to satisfy all equations in ([T6]) except those with n = 0, 1, . . . , A/", iV + 1. Applying 
([15]) to ([15]) also shows that formulae JT5J) are applicable at n = 1,7V. Homogeneous Dirichlet boundary conditions 
flno-i = &n -i = a ^o+i = ^m +i = ( a g ain for simplicity/specificity) imply that v - uj = a no+ i/a no = a mo _i/a mo 
and — (y + cj) = ^n +i/^n = ^m -i/^m which when combined with the ansatz ([T8]) gives 5 = —Ae lK ( 2n °~ 2 \ 
D = _Ce^( 2mo+2 ), P 7 = -^V^ 2 ™ " 2 ), D' = -CV^ 2mo+2 ). Thus we have 

( Ac n n < 1 f n < 1 , , 

an ~ \ Cd n n > N ' °" ~ \ C 7 < n> N [ZU) 

where c n is defined in terms of z = e m as before and c^ = (z') 77, - (V) 2no-2-n , ^ = _ ^2m +2-n £ or z / _ e */^ 
When combined with the dispersion relations ([T9]) . the condition ([20]) is enough to satisfy all except the 27V 
equations in ([T6]) associated with n = 1, . . . , N. Again we separate the remaining parts of the computation into two 
sections. 



1. Oligomer of length two 

For N = 2 the remaining four equations in ([TB]) for n = 1, 2 can be written as Pv = for P 

Pi = 



Pi Pi 
p 3 Pa 



and 
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(21) 



a 2 (^Y 

is a vector of length eight with each of a n , 6* restricted to the index < n < 3. By ([20]) . v can be 



v — 
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FIG. 4: The left set of panels corresponds to the VT- symmetric dimer with N = 2, ai,2 = 1 and Vi — ij — —Vi while 
the right set of panels to the case of the trimer with TV = 3, 0:1,2,3 = 1 and V\ = i*y = —Vz,Vi = 0. Each set contains a 
(top) plot of linear stability eigenfrequencies of the stationary solution corresponding to T = 0.7 and ko = 2.5 as computed 
on a lattice of length 20. Plots of the eigenvalue v indicating agreement between numeric (circles) and exact (x's) results 
are shown for the dimer with 7 = .5 (bottom left) and 7 = 1.5 (bottom right), and for the trimer with 7 = .5 (bottom left) 
and 7 = 1.5 (bottom middle) and 7 = 1.75 (bottom right). For the trimer, when 7 is small the dominant pair of (unstable) 
eigenfrequencies of negative imaginary part increases in magnitude as 7 increases until 7 ~ 1.15. At this point the complex 
pair recedes and a single dominant purely imaginary emerges at 7 ~ 1.3. This single eigenfrequency increases in magnitude 
until 7 « 1.49 when it begins to decrease. At 7 ~ 1.62, a new dominant complex pair increases in magnitude and continues 
to increase as 7 increases. 
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FIG. 5: Profiles of the moduli of unstable eigenvectors |a n |, |6*| are shown for the nonlinear dimer. The left four plots 
show eigenvectors for 7 = .5 and the right four plots for 7 = 1.5. These eigenvectors correspond to eigenvalues with negative 
imaginary part which are seen in Fig. [4] Again, numerical computations are shown to agree with the exact (x's) results and 
the length of the lattice is 20. Notice the localization in both components of the eigenvector in the case of the VT- symmetry 
breaking case of 7 > 1. 



The solvability condition is then the pair of equations 

det(PQ) = 0, z+^+u = -^z'+^+J S j (23) 

where the second was obtained from the dispersion relations ([T9]) . Solutions of ([23]) where either of z or z' is ±1 can 
be disregarded since the resulting system Pv = has no non-trivial solutions. The remaining solutions appear in 
quadruple sets (z, ^7), (-,z f ), (-, jr) so that the total number of relevant solutions of ([23]) is four times the 

length of the lattice. Thus by the calculation yields a number (counted with algebraic multiplicity) of eigenvalues 
v equal to the length of the lattice. 
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FIG. 6: Profiles of the moduli of unstable eigenvectors \a n \, \b^\ are shown for the nonlinear trimer with lattice length 20. The 
left four plots show eigenvectors for 7 = .5 and the right six plots for 7 = 1.5. These eigenvectors correspond to eigenvalues 
with negative imaginary part which are seen in Fig. 2] The agreement is similar to those of the earlier figures. 



2. Oligomer of length three or greater 

Similarly to the strategy for the linear case, for TV > 2 we begin by examining the 2 (TV — 2) equations in (p"6|) 
with n = 2, ...,7V — 1. Using ([20]) to rewrite the four quantities ai, <2at, ^i, fr^ one obtains a non- homogeneous 
square system where the 2(N — 2) variables a2, . . . , ajv-i, &2> ■ ■ • > ^jv-i can ^ e computed as linear combinations of 
A, C, A! , C' , under appropriate conditions of invertibility of the coefficient matrix. 

In the case of N = 3 the 2 (TV — 2) = 2 equations can be written as 



1 0J-V-V2 -2a 2 \^ 2 \ 2 1 
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(24) 



where the column vector is length six with 1 < n < 3. We move ai, 03, &*, 63 to the right-hand-side of and use 
(1201) to obtain 
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Of course, ([25]) is the nonlinear analogue of (|9j) and the values of z which are found by the remaining parts of the 
computation are such that the inverse matrix in ([25]) exists. The four remaining equations in ([T6]) with n = 1, 3 can 
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is length ten with < n < 4. By ([20]) and ([25]) we write 
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FIG. 7: The left set of panels correspond to the T^T-symmetric dimer with N = 2, 01,2 = 1 and Vi = ij = —Vi while the 
right set of panels to the case of the trimer with N — 3, 0:1,2,3 = 1 and Vi = 27 = — V3, Vi = 0. Each set contains contour 
plots of extremal stability eigenfrequencies \Im(is)min\ for extended solutions as in (|28[) on a lattice of length 20 with 7 = .5 
(top left), 7=1 (top right), 7 = 1.5 (bottom left) and 7 = 2 (bottom right). 




FIG. 8: The left set of panels correspond to the Hamiltonian dimer with N = 2, 01,2 = 1 and V\ = 1 + 8, Vi = 1 — S while 
the right set of panels to the case of the trimer with TV = 3, 0:1,2,3 = 1 and Vi = 1 + 5, Vi = 0, V3 = 1 — S. Each set contains 
contour plots of extremal stability eigenfrequencies \Im(u)min\ for extended solutions as in ([28]) on a lattice of length 20 with 
8 = (top left), 5 = 0.5 (top right), 5 = 0.75 (bottom left) and 5 = 1 (bottom right). 



The solvability condition is the same as ([23]) but with P and Q as defined in ([26]) and ([27]) . respectively. Again, 
solutions with either z or z' equal to ±1 are not relevant. 



B. Numerical results 



As mentioned in the Introduction, one of the motivations of the present work is to study the stability of scattering 
solutions. In particular, we focus here on of the class of plane wave solutions of ([T3]) of the form 



( R ikon + Re ~ikon n < x 
^ = { Te ik n n > 7V (28) 

with Ro,R,T G C representing the incident, reflected and transmitted amplitudes, respectively and ko > is the 
wavenumber such that uo = —2 co s (fen V. For each given values of fen,T one can compute ifj n by repeated application 
of the backwards transfer map 0, [ToMl2j 

^n-i = -^n+i + (V n - uj + a n \^ n \ 2 )^ n (29) 

which is a rearrangement of ([T3]) . For short oligomers, like the ones we are dealing with here, the ^ n can be evaluated 
analytically [HI, HiJ . 
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Fig. [4] shows an example of the the eigenvalues v computed with the method discussed in the previous section. 
As a check, we also report eigenvalues computed by numerical diagonalization of the matrix F in (fT6|) evaluated by 
substituting the computed ijj n in ([T7|) . The corresponding eigenvectors for the case of the dimer and the trimer are 
shown, respectively, in Figs. [5] and [6] 

It is clear that past the critical point of the VT phase transition, the relevant eigenvectors become highly localized 
(contrary to what is the case with the more spatially extended eigenvectors before the transition). These eigenvectors 
are responsible for the rapid growth of the norm density at the gain node observed previously [24] . The analysis 
also confirms the finding from our previous work [24[ that higher values of 7 and T correspond to more unstable 
solutions. We partially capture this phenomenon in Fig. [7J which contains a systematic two-parameter diagram 
of the dependence of the growth rate of the corresponding most unstable eigenstates. In addition to the stronger 
instability for higher T observed in the figure (and also for higher 7 in Fig. [4]), we observe the asymmetry of the 
relevant growth rate for transmission in the direction of positive vs. negative fco (as in pjl [24| we have adopted 
the convention that — fco labels a right-incoming solution with wavenumber fco). Among these two directions, it is 
intuitively clear and numerically confirmed that larger growth rates are incurred for waves that first encounter the 
gain site. For a fixed 7, the transition from right-incoming (—fco < 0) to left-incoming (fco > 0) waves may be 
achieved computationally by allowing the wave number to stay positive and instead flipping the potential V n from 
left to right. For the VT case this amounts to a continuous change in fco and a discontinuous change in the fixed 7 
value from negative to positive. This explains the discontinuity apparent in Fig. 

Finally, we touch upon the Hamiltonian case namely the case of real- valued V n for the nonlinear problem [l9| [once 
again we use a dimer with V\ = Vb(l + 5), V2 = Vb(l — 5) and a trimer with V\ = Vb(l + 5), V2 = 0, V3 = Vb(l — 5)]. 
This is captured in Fig. [8] In the Hamiltonian case the instability for the dimer associated with the particular values 
fco = 2.5 and T = 0.7, for example, increases as S increases from zero to S ~ 1.35 then the instability decreases for 
increasing 5. Additional oscillations appear in Figure 8 as an artifact of the short lattice length which we used to 
obtain the higher resolution with respect to the axes of Figure 8. These extra fluctuations disappear for the dimer 
as the length of the lattice becomes large. Nevertheless, the persistent characteristic corroborating the earlier work 
of [l9|| identified here concerns the asymmetry of transmission between incoming from the left (fco > 0) and from the 
right (fc < 0). 

IV. CONCLUSIONS & FUTURE CHALLENGES 

In this work, we have presented a methodology for addressing the spectral analysis of linear (or linearized) chains 
that possess an embedded oligomer defect. The technique, which draws a series of direct parallels with the approach 
used in [6|, consists of the solution of the linear problem (with its boundary conditions) on the lattice in which the 
defect is embedded. Subsequently these two linear solutions (on the left and right of the defect) are used as boundary 
conditions of the embedded region. Nontrivial solutions require that the determinant of a suitably defined matrix 
vanishes, ultimately leading to a polynomial problem for the eigenvalue, of the form given by ([23]) . Based on the 
solution of this resulting polynomial (which can be written fairly explicitly e.g. in the dimer case, see Eq. [7]), we can 
obtain the eigenvalues and reconstruct the eigenvectors of the full problem. 

We have applied the approach to genuinely linear problems with embedded complex oligomers, and, in the sec- 
ond part of the paper, to solve the stability problem as it arises from linearization around (analytically known or 
numerically computed) extended solutions with an embedded nonlinear defect. Finally, we calculated the spectrum 
and eigenstates in a number of cases of interest for T^Tsymmetric and Hamiltonian potentials. Our semi-analytical 
approach shows excellent agreement with numeric computations of the spectrum. 

The instabilities found in the nonlinear cases are generically oscillatory, with unstable eigenvalues having a non 
zero real and imaginary part and thus leading to oscillatory growth. The question, however, of the dynamical 
development of the instability and of the ultimate fate of such solutions is still in many respects open to more 
systematic studies. In the Hamiltonian case, there is numerical evidence that the instability leads to self-trapping 
of energy at the oligomer [25]. On the other hand, the dynamics of the VT case is fundamentally different in the 
instability development, usually ending in an indefinite growth at the gain site [24j . 

An especially interesting direction would be to consider generalizations of the present problem to higher dimensional 
settings. A priori, the solution of the linear problem is available in this case as well, e.g. for a square defect. Yet 
the matching of the solutions from the four relevant directions in the two-dimensional case presents considerable 
challenges. Generalizing the approach for both linear and nonlinear problems in such a case would be a particularly 
interesting topic for future studies. Another interesting issue to examine, even in the one dimensional setting, is the 
effect of boundary size on the eigenvalues of the problem (especially, so the PT-symmetric one). Both our earlier 
work [24] and that of others [26| have suggested interesting effects stemming from the finite size of the domain (and 
the associated boundary conditions) within which the embedded defect lies. Studying such effects and their scaling 
over the domain size would be another interesting direction for future work. Relevant themes are currently under 



consideration and will be reported in future publications. 
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